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Abstract. In a recent contribution, Dotsenko establishes a Fredholm determinant formula 
for the two-point distribution of the KPZ equation in the long time limit and starting 
from narrow wedge initial conditions. We establish that his expression is identical to the 
Fredholm determinant resulting from the Airy2 process. 

1 Introduction 

The Kardar-Parisi-Zhang (KPZ) equation is a stochastic evolution equation for a height 
function h(x,t), which models a surface growing through random deposition [1]. It was 
noted early on that the iV-th exponential moment of the height function can be expressed 
through the propagator of N quantum particles interacting with a short range attractive 
potential [2]. In one spatial dimension, this would be the attractive Lieb-Liniger Hamilto- 
nian. For its propagator we have now rather concise expressions, which allows one to carry 
out sophisticated replica computations so to obtain exact distribution functions for the 
height [3-11]. The various cases are distinguished by the initial conditions (narrow wedge, 
flat, two-sided Brownian motion in space variable x), by the location of the space-time 
reference points, and by possibly taking a scaling limit (e.g. t — > oo). 

In the recent contribution [12], Dotsenko develops a novel summation formula, which he 
uses to study narrow wedge initial conditions. In [13] he considers jointly h(0,t),h(0,2t), 
while in [12] he considers jointly h(x\, t), h(x2, t), both in the scaling limit t — > oo. It is 
this latter case which is of interest to us. Dotsenko obtains a Fredholm determinant for 
the joint distribution. Based on universality, Dotsenko's expression should agree with the 
corresponding Fredholm determinant of the Airy 2 process, but merely a direct inspection 
does not suffice to confirm the conjecture. 
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In our note we close the gap and prove the desired identity. As a by-product, the 
factorization assumption in [10] can be understood a posteriori as omitting the "cross 
terms" of Dotsenko's summation formula. 



2 One-point distribution 

The notations in [10,12] differ from each other. To make our contribution self-contained and 
to fix the notation, we briefly recall the replica computation for the one-point distribution. 
In dimensionless form the KPZ equation reads 

d, ifdhV ld 2 h 

dt h =2{d- X ) + 2dx~ 2+ ^ (2 ' 1} 

where h = h(x, t) is the height function at location x G R and time t > and rj is 
normalized Gaussian white noise. Through the Cole-Hopf transformation, Z = e h , (2.1) 
turns into the stochastic heat equation, 

We will consider only the narrow wedge initial condition Z(x, 0) = 8(x). 
Let us set 

h{x,t) = -±t + lt h{x,t), (2.3) 
7 t = (t/2) 1 / 3 , and introduce the generating function 

G ttX (s) = ( exp ( - e^*'**- 8 ) ) = < exp ( - e ( ' /24) - s Z(x, *)) ) . (2.4) 

Since lim^oo exp(— e~ ltX ) = 0(x), 9 the step function, 6{x) = 1 for x > 0, 9{x) = for 
x < 0, it holds 

lim G t2 ^J lt (s - x 2 )) = lim (6(s - h(2rfx,t) - x 2 )) = lim Prob\h(2-f 2 x,t) + x 2 < s] . 

(2.5) 

Gt, x ( s ) is a moment generating function for Z(x,t), 

oo / —s\N 

G t M =zZ- z ^( z ^ N > tN/2A - ( 2 - 6 ) 



N=0 



The N-th moment, (Z(x,t) N ), can be written using the eigenvalues and eigenfunctions of 
the 5-Bose gas with iV particles, 

(Z(x,t) N ) = <x|e-^*|0> = ^{x\^ z )^f g \0)e- E '\ (2.7) 
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where I X ^ — | X . • • • . X ), iV-times, H N is the Hamiltonian of the £-Bose gas of N particles, 
and z labels the eigenvalues E z and eigenstates Using the usual parametrization of 

the eigenstates [3,7], for fixed N one writes 



N , M , M 



2^ J 2^ a 24 



{x\* z )($ z \0) = |(0|^)|V^ W , (2.9) 



|(0|^)| 2 = iV!det( ^ -] , (2.10) 

2 ("a + np) + \{q a - qp) 



M 



TV . M oo 



M i 11 / 2 7T S.=i"., 

M=l 0=1^=1^°° 



Hence 



G tie ( s ) = ( -~ 1 - > JJ / ^ ^_ 1 ^n a -l e n3t/24+n a (i ?a:E -| g 2 t _ s) 

M=0 M! a=l 2?r n a =l 

M 



x det y- (2.12) 

z{n a + rip) + i{q a - qp) J a 



with the understanding that the term with M = equals 1. 
Using 

1 



\{n a + np) + i(q a - qp) 
and a simple identity for determinants, one has 

M 



i>oo 

/ duj a e~^ (na+n ^ +i(qa - q ^ a (2.13) 
Jo 



-»*/(»)) = E^-n/ ^/^Et-^r- 1 

M=0 a=l"' u ^ 00 n a =l 

xe «i*/ 24 +™c«(igc t 2:-5i3a*-*+ ;c2 / 2 *) (Jet | < e -5n Q (wc«+w (3 )-iij ct (w ct -w (9 ) N \ M ^ -j^N 

V / a,/3=l ' 



We rescale s — >■ 7tS,a; — >■ 2rfx,ou a — >■ ^t^a^Qa —> Qa/lt- Shifting g a — )■ g a + ix and 
noting that e^" - ^ inserted in the determinant does not change its value, we arrive at 

G t ^Ms - x 2 )) = det(l - , (2.15) 
where the kernel of K t s is given by 

^ e i ? (u ) '-u,)^ ( _ l) n-l e n3 t /24- 7t n(^ + i(a,W)- S )_ (216) 
2ff n=l 
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Using 



/oo 
dyAi(y)e Xny (2.17) 
■oo 

with A = t 1 / 3 /2 = 7t/2 2//3 and shifting y -> y + g 2 + |(u; + u'), one obtains 

K t ^,u')=2^J™ ge^'-) | dyAi(2 2 / 3 (y + ? 2 + ±(u; + u/)) 



oo 

< 

n=l 



^(-l)"-^-^"^) . (2.18) 

Summing over n and using the identity 

2 2/ 3 r ^ e 2iw Ai(2 2 / 3 (g 2 + a;)) = Ai(x + y)Ai(x - y) , (2.19) 

one obtains 

/°° e 7t(j/-«) 
dj/Ai(o; + y)Ai(u/ + y) - - . (2.20) 

In the t — > oo limit, 

P 7t(y-s) 

% - s) (2.21) 



1 + e^~ s 

and the kernel of K tiS turns into the kernel for the GUE Tracy- Widom distribution. 

3 Two-point distribution and factorization assump- 
tion 

Next we study the joint distribution of h(xi,t), h(x 2 ,t) for x 2 — X\ = 0(t 2 ^ 3 ) in the limit 
t — > oo. Following the strategy of Section 2, we consider the generating function 

G t , xl , X2 ( Sl , s 2 ) = <exp ( - e*** 8 *'')-! - e^^'*)- 8 )) 

= <exp ( - e^- Sl Z( Xl ,t) - e^- s *Z(x 2 ,t))). (3.1) 

It holds 

Inn G ti2j 2 Xu2l 2 X2 (^ t ( Sl - x 2 ),7i(s 2 - x 2 2 )) 

= lim Prob [^(xi, *) + x\ < s 1} h(x 2 , t) + x\ < s 2 ] . (3.2) 

By expanding in e~ Sl , e~ S2 , 

G t , xl , X2 s 2 ) = jrjt t?pLt£pl (z( Xl , t) L Z(x 2 , t) R )e^ R ^ , (3.3) 

L=0 R=0 
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where the term with L = = R equals 1. Expanding in the eigenstates of the 5-Bose gas 
one arrives at 

( Z ( Xl ,trz( X2 ,t)«) = ^i<°i^>i 2e " g ' t (gl, 7oitf^ > ' (3 - 4) 

where \x±, x 2 \ L,R) = \x±, . . . , xi, x 2 , ■ ■ ■ , x 2 ) with L times the argument xi and R times 
the argument x 2 . 

By Dotsenko's summation formula, the last factor in the summation can be written as 
(x 1 ,x 2 ;L,R\V z ) _ fL + RY 1 



(0\% 



m a +l a >0 ce=l 

.^^.^.-^^^^^^^^^^ x CTf (3 . 5) 

where x = x 2 — X\ and the summation over m a , l a is for all nonnegative values except for 
m a = = l a . The factor CT denotes the cross term which involves a product over all pairs 
a,/3, explicitly given by (cf. (26) in [12]) 

CT = (3.6) 
tt £[1 + \{m a + mp + l a + lp) + i(q a - qp)]T[l + \{-m a + m p + l a - Zg) + i(q a - q )} 

a-Tp H 1 + |( _m « + m /3 + l <* + W + Kfe - g/9)]r[l + |(m Q + 772,3 + / a - ^) + l(<7a - ' 

Substituting in (3.3), we obtain 

M=0 ^l" 7 - 00 m a +l a >0 V 7 

/ 1 \ M 

x det 



|(m a + l a ) + \(m p + l p ) + i(g a - qp) J a ^ =1 



xe ( m a+i a ) 3 t/24-|(m C( +« a )g^_ maS1 _; QS2 +i xirnaga +i X2 ; aga _ mc( ; C(:r /2 x ^ 

On the other hand, the expression (4.8) in [10] after factorization reads, upon replacing 
x — > x/2 there to match with our notation 1 , 



^tfxi .xo ( S l 



°° ( ~i\M roc i 00 / I 

.*) = E^-n/ ^EHr'^ L + i n / +i(g - qb 

M=0 m - a=l J-oo ^ na=l \ 2 n a + 2 np -i-H9a %> / a ,p=\ 

M 

X IT e"°' , ' 24 ~5 n ° , ° i e"^ 81 * ^gi^ifc-si _|_ gia;2(?c«-S2^ ri « ^ ^g g^ 



^ote that the definitions of ir in [12] and [10] differ by factor of 2. Also 1/2 in (3.21), (3.23), (3.27) 
of [10] should rather read 1/4. 



where di = d Sv i = 1,2. Binomially expanding the last factor on the RHS, we note that 

_ \ ^ e ~\xm a (n a -m a ) I n a \ e m a (ixiq a -si) e (n a -m cl )(ix2qa-S2) (3 9) 

^0 " W 



Rewriting the summation over n a in (3.8) by setting l a = n a — m a , this expression is the 
same as (3.7) except for the last factor CT. Hence the factorization assumption in [10] 
amounts to set the cross term CT = 1. 



4 Dotsenko's formula yields the Airy 2 process 

In [10] it is established that under factorization the limit t — > oo of the expression in the last 
section yields indeed the two-point distribution of the Airy 2 process. On the other hand 
in [12] Dotsenko argues that, after a suitable analytic continuation, CT — > 1 as t — > oo and 
then writes a Fredholm determinant expression for the limiting two-point distribution. One 
expects that Dotsenko's formula should coincide with the two-point function of the Airy 2 
process. In this section we establish rather straightforwardly the equivalence between the 
expression (51) in [12] and the t — > oo limit of (4.19) in [10]. 

Substituting s { 7^, x { -> 2rfx i: % — 1,2, cf. [10] Eq. (2.31), (4.19) in [10] turns into 

G ?2 7t V,2 7 f, 2 (^i - 7*(*2 - 4)) = det(l - AO , (4.1) 
where the kernel of N is given by 

/oo 
due -2 Xdl d2 e ^ldl+^d 2 Kd zl -d Z2 )^ t{u _ Si)>7t(u _ Sa) ) 
-00 

xAi(zi +M)Ai(^ 2 +u) • (4.2) 

Here 

$(U,V)= eU + eV = Y (-iy+™-lfl + ™\ mu+lv 

l+m>0 V 7 

Therefore 

lim $(7 t (u - si),7 t (u - s 2 )) 

t^oo 

= 0(u - si) + 0(u - s 2 ) - 0(u - si)0(u - s 2 ) =: $oo(si, s 2 , u) (4.4) 

and hence 

lim Prob[^(27 t 2 xi,t) + < s u h(2rfx 2 ,t) + x\ < s 2 ] = det(l - iV^) , (4.5) 

t— >OD 
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where 



/oo 
^ e -^ a2 e-^ ai+3:2% )^-^^ 0O ( Sl , S2 ^)Ai(z 1 +«)Ai(z 2 + W ). 
-oo 

(4-6) 

The derivative with respect to si, s 2 m &y seem rather formal, since it is acting on a step 
function. However, in later sections of [10] there are equivalent expressions for N^, for 
which such difficulties do not arise, compare with Section 5. 

Turning to Eq. (51) in [12], with the understanding that one should replace x — > 2 4 ^ 3 x, 
s — > 2 2 / 3 s to conform with the notations up to now, the probability in (4.5) is written as 
a Fredholm determinant with operator kernel 

A(u 1 ,u 2 ) 

C 00 f°° da 

= dy ^-Ai(y + q 2 + s 1 +u 1 +u 2 -iqx 1 )e iq{uJ2 - UJl) 
Jo J -oo 2tt 

+ r d V ? Ai {y + q 2 + s 2 + co 1 +u 2 + iqx 2 ) e 1 ^^ 

Jo J-oo 21T 

- f dy %L (27r)- 3 / 2 | d^d&d&Aiiy + q 2 + Wl + u 2 - i^)^ 2 ^ 

xe ~h^~^/^) 2 -^-^/^) 2 -me(y - si + ZiV^)8(y -S2 + &Vx) 

= J ^ J dyAi(y + q 2 + u 1 +u 2 )e iq{uJ2 - UJl) {d(y-s 1 + iqx 1 )+9(y-s 2 + iqx 2 ) 

-(2tt)- 3 / 2 J ^irf6^3e"^ {f?+522+€|) e- (5l+if3) ^ 1 e-« 2+i ^)^ 2 

x 9 (y - si + iqxi) 6(y- s 2 + iqx 2 ) } . (4.7) 

Here the step functions with complex arguments are to be understood as resulting from a 
limit in analogy to (2.21). The last term can be rewritten as 

(2vrr 3 / 2 j d^^^e-^^^-^^^e-^+^^eiy -s 1 + iq Xl ) 

x9(y-s 2 + \qx 2 ) = e- xdld2 9(y - Sl + iqx 1 )9(y - s 2 + iqx 2 ) . (4.8) 
Proof. It is easy to see that 



g— \xdid2^— msi — ls2 g— msi— ls2 — \xral 



= (27r)~ 3/2 J ^ 1 ^ 2 ^ 3 e~5(??+«2+?3)+ m v^6+'v^6+«(m+0v^6 e -msi--«s2 

= (2n)- 3 / 2 J ^ 1 ^ 2 ^3e-|(^+^+^) e -^+^)Viai e -«2+^3)^ e -m S1 -i S2 _ ( 4 9 ) 

Multiplying both sides by (-iy+rn-i e mu+iv+i(mx 1+ ix 2 ) q (m+ij ? substituting Si -> 7^, x t ->■ 
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and summing over l,m, we obtain 
e- xdld2 $(7 t (u - si + ixig),7t(« ~ s 2 + iar 2 g)) 

xe -(6+i6)v / 5a le -(6+i?3)^ 2$ ^ 7 ^ M _ Si + ixi? ) )7t ( M _ g 2 + ix2 g)) . (4.10) 

Taking the limit t — > 00, we arrive at the desired equality. □ 
Using Eq. (4.8) one writes 

A( Wl) w 2 ) = J°° dyAi(y + q 2 + u 1 + uj 2 )^~^ 

xe -l8iai $ M (si - Lrig, s 2 - i^ 2 g, 2/) 
= 2 ^°° g e i^/3((. 2 -. l) -K9 1+ ^ 2 )) j d y k i(2 2 '\y + q 2 + 2- 2 l\uj l + uj 2 )) 

xe- x9 ^ 00 (s u s 2 ,2 2 / 3 y). (4.11) 
We switch to x, s as in (4.6) and use (2.19) to obtain 

A(u u u 2 ) = 2 1/3 J dyk\(y + 2 1/: V - {x^ + :r 2 <9 2 )) Ai(y + 2 1/3 w 2 + (x^ + x 2 d 2 )) 

xe-^Q^s^y). (4.12) 

Since $oo can be written as a limit of a sum of exponentials, see (4. 3), (4. 4), the derivatives 
<9j, i — 1, 2, inside the Airy functions can be regarded as numbers. Finally one can rewrite 
(4.12) as 

A{u u u 2 ) = 2 1 /3y^ e -^ 1 a 2e - 2 - 1 /« (2;i a 1 +x 2 9 2 )(^ 1 -^ 2 ) 

x Ai(y + 2 1 /V)Ai(y + 2 1 / 3 W2 )$ 00 ( Sl , s 2 , y) , (4.13) 

which is the same as the kernel A^ of (4.6) with the identification Zi = 2 1 ^ 3 u i ,i = 1,2. 

We draw one conclusion from our computation. Both in [12] and [10] one has to linearize 
the factor e~ xml ) see (4.9). [12] uses a threefold Gaussian integration, whereas [10] uses the 
differential operator of (4.9). If already in [12] Eq. (34) one would use the latter method, 
then A^oo would appear without further efforts. 



5 A simplified proof for the Fredholm determinant of 
N to determine the Airy2 process 

In [10] it is proved that the Fredholm determinant of A" equals the one of the two-point 
distribution of the Airy 2 process as t — > 00. The arguments of [10] are somewhat involved. 
Here we provide a more straightforward proof. 
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We start from an expression Gf x X2 (si — x\/2t, s 2 — x 2 /2t) = det(l — N), which follows 
from multiplying out the last factor of (3.8), with the kernel of N given by 

oo n , x 

N(z,z') = l { ^ >0} ^^(-l)MHe-( Sl -^/ 2 ') fc -( S2 -^/ 2 ')(«- fe ) 

71=0 fc=0 ^ ' 

xe-i x * M) J duki(u + z + y)M(u + z' - y)e^ tUn (5.1) 
with y = (xik + a; 2 (n — k))/2^ t . Let us rewrite the integral over u as 
J dudv5(u - v)e ytUn Ai(u + z + y)M(v + z - y) 

= J dudv5(u - v)e^ nw e^ n ~ k) Ai(u + z + (n - k)x/2-f t )e ,tvk Ai(v + z' + kx/2 lt ) , (5.2) 

where the integrals over u,v have both been shifted by w = {—x\k + X2{n — k))/2^ t - We 
notice the identity 

e a "Ai(w + A) = e - r(A+a2) e ar2 e- r ^e au Ai(t; + A - 2ra) (5.3) 

with the Airy operator H = — (d u ) 2 + u, see Appendix B of [10]. 
Proof. Since 

(H + A)Ai(u + A) = , (5.4) 

one arrives at 

e au Ai(« + A) = e^ H+x) e T ^ H+x) e au e~ r(H+x) Ai(u + A) . (5.5) 

Using 

u(t) = e TH ue~ TH = u + r 2 - 2rd u , 
(5.3) follows from 

e au Al(u + A) = e - T(x+a2) e aT2 e- TH e au e- 2Tad "Ai(u + A) . (5.7) 

□ 

With the help of identity (5.3) the kernel of N can be rewritten as 

oo n s \ 

N(z,z>) = l {2>2 , > o } ^^(-l)- 1 (^)e-- fc --^ fe )e^-'-^)/^ 2 

n=0 fc=0 ^ ' 

x J dudv{u\e- xH/2 ^\v)e lt{n - k)u+ltkv Ai{u + z)A\{v + z') . (5.8) 



Using the cyclicity of the determinant, det(l — AB) = det(l — BA), it holds that det(l — 
N) = det(l — L) with the kernel of L given by 



oo n / \ 

£(».«) = ED- 1 )- 1 UJ* 

n=0 fc=0 ^ ' 



k-s 2 (n-k) l u \ e ~ x 




xe 



Clt{n— k)u+jtkv 



/ Ai(w + z)Ai(t> + z)e 



xz/2-yf 



'Hz. 



(5.9) 



From this expression it is easy to see that the Fredholm determinant of L yields the one 
of the Airy 2 process for t — > oo, see the arguments in Section 2.4 of [10]. 

Acknowledgement. We are grateful to Victor Dotsenko for useful comments and for 
communicating his preprint [12]. 
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